VENUS ATMOSPHERE PROFILE FROM A MAXIMUM ENTROPY PRINCIPLE 
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The variational method with constraints recently developed by Verkley and Gerkema to describe 
maximum-entropy atmospheric profiles is generalized to ideal gases but with temperature-dependent 
specific heats. In so doing, an extended and non standard potential temperature is introduced that is 
well suited for tackling the problem under consideration. This new formalism is successfully applied 
to the atmosphere of Venus. Three well defined regions emerge in this atmosphere up to a height of 
100 km from the surface: the lowest one up to about 35 km is adiabatic, a transition layer located 
at the height of the cloud deck and finally a third region which is practically isothermal. 



1. Introduction 

Since the late 1950s, there have been hundreds of space 
exploration missions that have provided data on the com- 
position, structure and circulation of planetary atmo- 
spheres. In particular, the vertical pressure-temperature 
profile of the atmosphere of Venus was accurately ob- 
served in the Pioneer- Venus II mission in 1978 hi. In 
ref. 0, an averaged data of temperature, pressure and 
density of the Venus atmosphere as a function of altitude 
above the surface is given. In this paper we will use this 
information. 

Venus is the nearest planet to the Sun. The incom- 
ing solar energy, of the order of 2600 W/m 2 , is almost 
twice as great as on Earth. Venus possess an atmosphere, 
which is about a hundred times as massive as that of the 
Earth and its surface temperature reaches 730.fr. Car- 
bon dioxide is the major component (96.5 per cent), and 
3.5 per cent of N2 is the next most abundant species. 

In a recent paper by Verkley and Gerkema [3j, hence- 
forth denoted by VG, have proposed a first-principles 
variational method to generate a p — T profile for an 
atmospheric column. This profile is the result of maxi- 
mizing the entropy of the column, S, maintaining fixed 
its mass, M, enthalpy, H, and the integral of the po- 
tential temperature, L. This method is flexible enough 
to describe from an isothermal profile to an adiabatic 
one depending on the value of a unique parameter. VG 
applied these ideas to the Earth's troposphere and accu- 
rately reproduced the U.S. Standard Atmosphere. 

In their analysis, VG considered that the atmospheric 
gas is ideal, i.e. a gas that verifies the corresponding- 
equation of state, and besides, the specific heats, Cy ;P , 
are constant. When this second assumption is also ful- 



filled it is said that the gas is perfect. These two hy- 
potheses are quite reasonable when dealing with the at- 
mosphere of the Earth. For representative temperatures 
and pressures of the low Venus atmosphere, however, the 
temperature dependence of the Cy,p has to inescapably 
be taken into account. As our goal in this paper is the 
application of the VG ideas to the Venus atmosphere, we 
have had to extend their original variational method to 
take into account this dependence of the specific heats on 
temperature. Consequently we will first extend the VG 
method to ideal imperfect gases. As is well known, the 
temperature dependence of the Cy iP comes from the exci- 
tation of some internal vibrational levels of the molecule 
when the temperature is high enough. 

Throughout the paper we will consider that CO2 is the 
main chemical component of the Venus atmosphere. The 
chemical impurities coming from the N% produce negli- 
gible corrections in our analysis. Besides, we will not 
consider any correction coming from the interaction be- 
tween the molecules, that is, we will stick to the scenario 
of an ideal gas but with the afore-mentioned tempera- 
ture dependence of the Cy, P - Corrections of real gas in 
the Venus atmosphere have been considered but for 
our purposes here they are negligible. Moreover, we will 
consider a constant gravity for the entire column under 
study because its variation with altitude has no impact 
in the results. The verification that the equation of state 
of ideal gases is good to describe the gas of Venus can 
be done using the data of ref. [2J. There one can check, 
at each height, that using pressure and temperature as 
input, the value of density obtained from the equation of 
the ideal gases is, in all cases, equal to the experimental 
value with a discrepancy lower than 2 per cent. 

The index of the paper is as follows. In Section 2, we 
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will recall the meaning of the potential temperature, 8, 
for ideal perfect gases and its relation to the entropy. In 
Section 3, we will review the VG method and comment 
on the main results when it is applied to the Earth's tro- 
posphere. In Section 4, we will put in evidence the ne- 
cessity of considering an explicit dependence of the Cy, P 
with temperature for the CO2 molecule, to properly de- 
scribe the Venus atmosphere. In Section 5, we will cal- 
culate the equation of the potential temperature, i.e. the 
form of an adiabatic, when dealing with ideal-imperfect 
gases. In Section 6, we will extend the VG method to 
ideal-imperfect gases. In Section 7, this new method will 
be applied to the Venus atmosphere and the results will 
be commented on. Finally, in Section 8 we present our 
conclusions. 



2. The Potential Temperature and the Entropy for 
Ideal Perfect Gases 

The first principle of Thermodynamics says that 

6U = 5Q-pSV. (1) 

Throughout this paper, 8 indicates all the variations of 
thcrmodynamical quantities, whereas d will be used to 
refer to variations of quantities at different altitudes. 
For ideal gases, the equation of state, for 1 mol, is 



pV = R* T 



(2) 



(where R* = 8.3143 J mol -1 K^ 1 is the universal gas 
constant) and the specific internal energy (internal en- 
ergy per unit mass, noted with lower case letter) reads 



As it is standard, choosing the reference pressure p r 
1000 mb, the potential temperature 8 reads 



e = T 



1000 mb 
P 



(7) 



The potential temperature of a gas is therefore the tem- 
perature that it would take if we compress or expand it 
adiabatically to a reference pressure of 1000 mb |f| . Ob- 
viously, the choice p r = 1000 mb is only a question of con- 
venience when we are dealing with the low atmosphere 
of the Earth. In other words, the value of the reference 
pressure p r is arbitrary. 

The differential expression for the specific entropy of 
ideal gases is 



a r 5T 

OS = C-n — — 



(8) 



Because of the definition of 8 in Eq. J7J) , 8 In ( 
n8hip and consequently 



SlnT 



Ss = C p 8\n8, 



(9) 



and because we are considering that C p is constant 



s = C p ln0- 



constant. 



(10) 



Due to this direct relation, the potential temperature 8 is 
usually considered as the entropy of the gas in the physics 
of the atmosphere. It is apparent that in a process of 
constant 8, the entropy is also constant. 



3. The Verkley-Gerkema Method 



= Cy T + constant. 



(3) 



In an adiabatic process, 8Q = 0, using the first two equa- 
tions and remembering that C p = Cy + R* , we find 



r ST Sp 
C p —-R -. 



(4) 



Thus, two points (p, T) and (p r ,T r ) belonging to the 
same adiabatic trajectory fulfil 



Co 



ST' 



= R* 



Pr P' ' 



(5) 



In this case, the reference temperature T r becomes the so 
called potential temperature 8 corresponding to the ref- 
erence pressure p r . Note that in the previous equations, 
C p has been considered as a temperature independent 
constant because in this section we are dealing with an 
ideal-perfect gas. This last equation constitutes one of 
the forms of expressing Poissons equation for an adia- 
batic process. Explicitly it reads as 



Tp~ 



TrPr 



R* I C p . 



(6) 



VG have proposed a first-principles variational method 
to generate ap-T profile for an atmospheric layer. In 
this method, one considers a column of dry air with unit 
horizontal area in hydrostatic equilibrium, bounded by 
two fixed values of the pressure. At each altitude z, the 
local thermodynamics is related to p(z), T(z) and p(z) 
that can be inverted into z — z(p), T{p) and p(p), since 
dp/dz = —g p. The profile one is looking for comes out 
after the maximization of the entropy with the condition 
that the enthalpy and the vertical integral of the potential 
temperature, L, are kept fixed. As the top and bottom 
pressures, P2 and p\, of the column are specified, the total 
mass of air is also fixed. We denote by 
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the total entropy, mass, enthalpy, and vertical integral of 
the potential temperature of the air column, respec- 
tively. The requirement of L fixed is appropriate in the 
case of convective mixing. An analysis on this issue and 
its bibliographic antecedents are presented in ref. 

As said above, VG imposed the maximization of S, 
with M , H, and L fixed, namely 



5S + \5H + [i5L = 0, 



(12) 



A and /i being Lagrange multipliers. This means that one 
should maximizes the entropy with respect to the tem- 
perature profile of the atmosphere with the constraints 
imposed. The profile variation is consequently performed 
along the column through the variation of the tempera- 
ture keeping the pressure fixed at every altitude. 

From the maximization with constraints, VG deduced 



Introducing 



T \p 



= 0. 



(13) 



and denoting the temperature at po by To, the resulting 
p — T profile is 



T = T 



1 



1 + a {po/pY 



(14) 



This result is actually powerful; note that it is able 
to describe a wide range of profiles, from an isothermal 
profile, in the a — > limit, to an adiabatic profile in the 
limit of a — > oo. In this last case, po and To clearly be- 
come p r and 9 respectively. In the application of this 
result to a specific air layer, a is fixed by requiring that 
the ratio H/L, that derives from the profile, is equal to 
its real value in the air layer. Then, To is fixed by re- 
quiring that H (or L) is equal to its real value as well. 
Using this method VG obtained a profile for the tro- 
posphere that fitted remarkably well the U.S. Standard 
Atmosphere data. The same method has been applied to 
the Earth's mesosphere Q- 



4. Specific Heat of CO2 

The molecule of CO2 is linear and consists of a carbon 
atom that is doubly bonded to two oxygen atoms. Its 
molecular mass and corresponding specific gas constant 
are 




FIG. 1: Specific heat at constant pressure (C' p ) in units of 
R* as a function of temperature (T), of the CO2 molecule 
(continuous line) and its polynomial fit (dashed line). For 
details see Section 4 



type are expressed in the form: 



sinh 2 ( — — 
2T 



(15) 



The first term in this sum corresponds to translation and 
to rotational modes while the second term is due to vi- 
brational modes. The sum in the second term extends 
to the various characteristic temperatures, T, , of the 
vibrational modes. For CO2, these values are M 



T V (K) = 960, 960, 2000, 3380. 



(16) 



For practical uses, it is convenient to express Ea. l|15|l 
by a simple algebraic approximation. Thus, it is stan- 
dard to express the specific heat by a series of powers in 
T, with coefficients empirically adjusted. Three terms, 
namely 



C„ = A + BT + CT 2 , 



(17) 



are enough to describe the specific heat with an excel- 
lent accuracy in the range we are interested in. This is 
clearly shown in Fig.l, where we have plotted C p vs. T 
for COi. The continuous line corresponds to the "exact" 
values which have been obtained by plotting (|15J) with 
the insertion of the four characteristic temperatures of 
vibration as given in ljl6|l . The dotted line is the polyno- 
mial approximation (|17fl . the coefficients being 



A = 2.5223, 5 = 0.7710110" 



C 



-0.3981 10" 



M, 



co 2 



The units of A, B and C are R* , R*/K and R*/K 2 
U.Olgmor 1 ; R = M C o 2 *R* = 188.92 J kg^ 1 # respectively. From Fig.l it is apparent, first, that C p 

has a significant change in the range of temperatures of 
the Venus atmosphere, from 735 K at z — to 175 K 
at z — 100 km. Even if we restrict our analysis to the 



The various contributions to the specific heat, at con- 
stant pressure, of a gas constituted by molecules of this 
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lowest layer, i.e., from z — to z — 40 km, the interval of 
temperature is between 735 K and 385 K. In Fig.l one 
first observes that the change in C p in this interval of 
temperature is of the order of 50%, and secondly, that 
this polynomial fit is excellent. Thus, we will use H17|) as 
input to generalize the VG method and apply it to the 
atmosphere of Venus. 



5. Potential Temperature for a 
Temperature-Dependent Specific Heat 

The generalization of the VG method can start from 
Eq.Q) because it is applicable to any ideal gas. In this 
Section however, Eq.lJSJ for the adiabatic transforma- 
tions, adopts the form 

£ C P {T>) (A+BT'+CT' 2 ) ^ = R* J" *± 

(18) 

p r being the pressure level of reference and 9 the potential 
temperature characteristic of this adiabatic trajectory. 
We find 



P_ 

Pr 

(19) 

In Fig. 2, we have adjusted the data of the low atmo- 
sphere of Venus by means of this formula for an adiabatic. 
There is a generalized consensus supporting this behavior 
as can be read in ref . Jl| . It shows an excellent agreement 
along the lowest 30 km and a clear departure for upper 
heights. 

In the present case, the physical meaning of the po- 
tential temperature is the same as that of J2J), i.e. the 
temperature that the gas would take if we compressed or 
expanded it adiabatically to a pressure p r . The price to 
pay now is that given the coordinates (p, T) , the value 
of 9 has to be computed numerically by iterations. Of 
course, in the case B = C = 0, J7J and (|19() are the same 
equation. 

The differential expression for the specific entropy of 
ideal gases was written in (|5J) and now reads 





~B{T -9)~ 




~C(T 2 -9 2 )~ 




R* 


exp 


2R* 



Ss = C p (T) 6 -f^R* S f 



(20) 



For practical purposes, let us define a new magnitude 
r, with physical dimensions of temperature, in the fol- 
lowing way: 



with 



C° = C P (T ) ; r(T ) = T . 



(21) 



(22) 
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FIG. 2: Pressure-temperature profile of the atmosphere of 
Venus. The data (black dots) have been taken from Table 1, 
and the fitting line is an adiabatic as given by Eg. 11911 . 



In fact, this new variable r,that reduces to the stan- 
dard temperature as soon as C p is constant, allows one 
to treat this problem in exactly the same way as in the 
case of the ideal perfect gas. 

After specifying C P (T) in a form like 117(1 . and by the 
integration of l|21l) we obtain 



St 

T 



{A 



ST 



A 



In 



To) 



BT + CT 2 ) 
B 
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T 

To) 



T - Td 



and the result is 
\T CI 



T 

7f exp 



f(T-To) 



This is the relation between r and T 
fulfils the two agreed upon conditions: t(To) 
t(T) = T for B = C = 

Using (ED in JU) 



C_ 
2A 

(24) 

This formula 
Tq and 



5s 



St 



R* Sp 



(25) 
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and in a similar way to what was done to arrive to (I10|) . 
we obtain 



s = Cp lnf + constant, 



with 



and 



r = r | — 



R* 



(26) 



(27) 



(28) 



It is now clear from this result that in an isentropic pro- 
cess, f is also conserved. The use of lf2*ri)) will permit us 
to extend the VG method quite easily, at least from a 
conceptual point of view. 



6. Generalization of the VG Method 

The formal analogy between these new magnitudes and 
those appearing in Section 2, f «-> 9, k° «-> k , permits 
us to extend the VG method reviewed in Section 3 in a 
quite natural way. The entropy, the enthalpy and the L 
function of an atmospheric column per unit section are: 



S< = _p / lnf dp, 

9 Jp 2 



(29) 



1 



H = - I hdp=- CJT)Tdp, (30) 
9 J P2 9 J P2 
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FIG. 3: Pressure-temperature profile of the atmosphere of 
Venus. The data (black dots) come from Table 1 and the 
theoretical results (continuous line) have been obtained using 
the extended VG method. 



L = — / fdp. 

9 J P 2 



(31) 



The condition of maximum S, with H and L constants, 
is expressed as the extreme of the functional ^ under the 
variation of the profile T(p), in the following way 



5* 



pi fiO 

p 



rPi i rPi 

5 [lnf + /i f ] dp+X - S[h]dp = 0. 
1 pi 9 J pi 9 J pi 

(32) 

To go into the details of this variational problem, re- 
member that 

t = t(T) ; h = h(T) ; f = f(r,p), (33) 

Then, the condition 8^ — (p fixed) implies that 



1 \ 6f X Sh 

f + n 5f + C°5f = °> 



and taking into account that 

5 lnf Sh 



TCl- 



(34) 



(35) 



one gets 



T 



+ X = 0. 



(36) 



Finally, the result of the extreme can be written as 



with 



T = T a 



X 



1 



1 + af/T' 



T ° = x- 



(37) 



(38) 



being formally identical to the VG result. Nevertheless, it 
does not provide an explicit expression of T as a function 
of p. This is due to the fact that f depends on r and 
through it on T, so that f/T is now a function of p and 
T. Consequently, the equation (|3Tjl has to be treated 
numerically. In fact, by using the method of Lagrange 
multipliers, the experimental data fix the constraints H 
and L. Then, the parameters To and a are determined 
in order to satisfy H and L. 
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7. Results of the Extended VG Method 

The performance of the extended VG method, sum- 
marized in Eq. l|37jl. when applied to the atmosphere of 
Venus is shown in Fig. 3. Note that the logarithmic scale 
used for pressure is more sensitive and hence puts more 
clearly in evidence the virtues or flaws of the fit. The 
characteristics of the Venus atmosphere allow us to well 
distinguish three regions in altitude. In each one of these 
regions we have applied the extended VG method by 
maximizing the entropy with both H and L constraints. 

In the lowest one, from z = to z ~ 34 km, we found 
that the value of a is 44 and the value of To = 735 K. 
This high value of the unique parameter of the model 
indicates the clear preponderant role of L with respect 
to H, as the dominant constraint in the maximization 
of S. From this fact we conclude that in this layer the 
profile is adiabatic. This clearly agrees with the result 
shown in Fig. 2, obtained by means of Ea. (|19|l . 

In the layer from z ~ 60 km to z ~ 100 km, the value 
obtained for a is 0.18 and To = 245 K This indicates that 
here H is the dominant constraint, and in consequence 
the profile of this layer is basically isothermal. 

The method assigns a — 3 with To = 444 K to the 
intermediate layer, a transition region, from z ~ 34 Km 
to z ~ 60 Km. 

8. Discussion and Conclusions 

Our purpose in this paper was to extend the first- 
principles VG method to atmospheric layers where the 
thermal deviation is so high that it is not reasonably 
possible to consider a constant unique value for the spe- 
cific heat in all the points of the layer. To implement 
this extension, we have had to define a new extended po- 
tential temperature, f. A potential temperature is, by 
definition, a generalized temperature that characterizes 
adiabatic processes, or in other words, that is conserved 
in processes where the entropy is constant. The rela- 
tion between the new f and entropy maintains the same 
standard form as that existing for ideal-perfect gases. 

^From the surface of Venus up to about 100 km, we 



have distinguished 3 layers. The extended VG method, 
for each layer, after fixing To to the corresponding initial 
temperature of the layer, deals with a unique parame- 
ter a to adjust the constraints. From the surface up to 
35 km, the method assigns a high value of a. This implies 
that the p — T profile here is adiabatic. There is no sur- 
prise in this result; an adiabatic lapse-rate was detected 
in the low atmosphere of Venus during the first obser- 
vations of the Russian Venera probes. This was verified 
later in other spatial missions [lj. In the layer from 35 
to 65 km in height a has an intermediate value that in- 
dicates a transition. Finally, from 65 to about 100 km, 
this method produces a low a, which is interpreted as 
an isothermal layer. This is the first time that the VG 
method identifies an adiabatic layer where it should; this 
concordance shows clearly the success of the method it- 
self and of the extension presented in this paper. 
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